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We report on the current status of knowledge on inclusive quarkonium production in high-energy photon- 
photon collisions. As a perspective for the LHC, we compute various production cross sections via direct photon- 
photon fusion in ultra-peripheral pp, pA and AA collisions at the LHC using the tree-level quarkonium amplitude 
generator MadOnia. 



1. Introduction 

In a few months the LHC will start operat- 
ing and provide physicists with a tremendous 
amount of new experimental information on the 
fundamental particles and their interactions. In 
this context, much is expected from the LHC 
to elucidate the still ill-understood quarkonium 
production mechanism (see |l|2|3j for recent re- 
views) . This is of particular importance in view 
of the possible suppression of quarkonium pro- 
duction after quark-gluon plasma formation in 
heavy- ion collisions. One of the questions which 
remains to be answered is in particular the rel- 
evance of color-octet transitions in different pro- 
duction regimes or equivalcntly of the 0(v) cor- 
rections in Non-Relativistic Quantum Chromody- 
namics (NRQCD) [4l5] , 

The large amount of quarkonia that will be 
produced in the pp, pA and AA modes at the 
LHC will allow for even more precise analyses 
than those performed recently at Run II of the 
Tevatron |6l7j . We expect cross section and po- 
larization measurements of inclusive, prompt and 
(hopefully) direct yields of the J/ip, ip', T, T' and 
T" mesons. In addition, analyses of their associ- 

* Research supported by the French Ministry for Higher 
Education and Research, CNRS-IN2P3 and ANR. 
tTalk given by J.P. Lansberg at the Workshop on "High- 
energy photon collisions at the LHC", April 22-25, 2008, 
CERN, Switzerland. 



ated production with a heavy-quark pair [819] or 
of other channels would certainly be as interesting 
and in some cases better suited to disentangle the 
color-singlet and color-octet transitions and to as- 
sess the importance of s-channel cuts |10lllj and 
of the rescattering with comovers [12113) . Re- 
cently, several quarkonium production cross sec- 
tions have been evaluated at next-to-leading or- 
der (NLO) of perturbative QCD |14|15|16j . com- 
plementing the pioneering work for color-singlet 
production in direct jp collisions [17j . The com- 
mon feature of these calculations is the signifi- 
cant size of the NLO corrections, in particular for 
large transverse momenta of the quarkonia. The 
color-singlet prediction can thus be brought con- 
siderably closer to the data. As of today, only 
the color-octet contributions to direct 77 colli- 
sions have been evaluated at NLO [14I15I I 1 ! 

Very little is known yet about NLO effects 
on the polarization of the produced quarkonia. 
However, the leading order (LO) NRQCD predic- 
tion including sizable color-octet contributions is 
known to be in conflict with the measurements 
e.g. in hadroproduction [19120] , A recent NLO 
evaluation for J/ip mesons produced directly as 
color-singlet states in gluon-gluon fusion changes 
the LO prediction from almost purely transverse 
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to rather longitudinal polarization at NLO [2T] . 
However, the corresponding unpolarized differen- 
tial cross section still falls short of the data, pre- 
venting any conclusive interpretation of this re- 
sult [16] . A forthcoming study of both the differ- 
ential cross section and the polarization at NLO 
complemented by dominant NNLO corrections 
might help to elucidate the situation further, at 
least for hadroproduction [22] , 

In the absence of an International Linear Col- 
lider (ILC), the LHC will, at least in the short 
term, present an almost unique environment for 
high-energy photon-photon collisions [23] ■ These 
tend to produce much cleaner events than those 
produced in inelastic hadron collisions. They may 
also help to validate the only experimental anal- 
ysis of quarkonium production in photon-photon 
collisions to date [24j . which clearly required im- 
portant color-octet contributions [25]. We there- 
fore focus here on quarkonium production in 
photon-photon collisions at the LHC. First, we 
recall the theoretical status along with more spe- 
cific results related to the LEP II analysis. We 
then present original results for ultra-peripheral 
pp, pA and AA collisions obtained with the tree- 
level quarkonium amplitude generator MadOnia 
[27] . before we present our conclusions and an 
outlook. 

2. Inclusive quarkonium production in 
photon-photon collisions 

Photons can interact either directly with 
the quarks participating in the hard-scattering 
process (direct photoproduction) or via their 
quark and gluon content (resolved photopro- 
duction). Thus, inclusive quarkonium produc- 
tion in photon-photon collisions receives contribu- 
tions from the direct, single-resolved, and double- 
resolved channels. All three contributions are for- 
mally of the same order in the perturbative ex- 
pansion and must be included. The J/ip mesons 
can be produced directly, or via radiative or 
hadronic decays of heavier charmonia, such as 
Xcj and ip' mesons, or via weak decays of B 
hadrons, where the latter can often be safely ne- 
glected. The cross sections of the four resid- 
ual indirect production channels may be approxi- 



mated by multiplying the direct-production cross 
sections of the respective intermediate charmo- 
nia with their decay branching fractions to J/ip 
mesons. 

Invoking the Equivalent Photon Approxima- 
tion and the factorization theorems of the QCD 
parton model and NRQCD, the differential cross 
section can be written as 

da(AB -> AB + H + X) = J dx A f l/A {x A ) 
x / dx B f 7/B (x B ) y~] dx a f a/l (x a ,M) 

J a,6,d 

x I dx b f bh (x b ,M)Y,(0 H [n}) 

J n 

x da(ab — > cc[n] + d), (1) 

where fj/A,B( x A,B) is the equivalent number of 
transverse photons radiated by the initial-state 
particles A and B, f a ,b/-i{x a ,b, M) are the parton 
densities of the photon, (O h [n]} are the operator 
matrix elements of the H meson, da(ab — > cc[n] + 
d) are the differential partonic cross sections, 
the integrals are over the longitudinal-momentum 
fractions of the emitted particles w.r.t. the emit- 
ting ones, and it is summed over a, b — J,g,q,q 
and d — g,q,q, with q — u,d,s. To leading or- 
der in v, we need to include the cc Fock states 
n = ^^M^^Pf if H = J/ip,ip> and 

n = 3 P^\ 3 S[ 8) if H = XcJ, where J = 0,1,2. 
With the definition / 7 / 7 (x 7 ,M) = S(l — x 7 ), 
Eq. ([1]) accommodates the direct, single-resolved, 
and double-resolved channels. The presence of 
parton d is to ensure that Pt can take finite val- 
ues. 

In Fig. [TJ the P^ distribution of e + e~ — * 
e + e~+J/ifj+X measured by DELPHI [24] is com- 
pared with LO NRQCD and color-singlet predic- 
tions 25J. The solid lines and shaded bands rep- 
resent the central results with MRST98LO par- 
ton densities and their uncertainties, respectively, 
while the dashed lines represent the results ob- 
tained with CTEQ5L parton densities. It is clear 
that at LO the DELPHI data favor the NRQCD 
prediction, while they significantly overshoot the 
CSM one. Note, however, that the associated pro- 
duction of J/ip mesons with additional cc pairs 
can be numerically important |26j . 
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NRQCD e+e ~~ > e+e x at LEP2 




Figure 1. The cross section da/dP^ of e + e — * 
e+e" + J/if) + X measured by DELPHI [24| as 
a function of P T is compared with the theoret- 
ical predictions of NRQCD and the CSM. The 
solid and dashed lines represent the central pre- 
dictions obtained with the ME sets referring to 
the MRST98LO (default) and CTEQ5L PDFs, 
respectively, while the shaded bands indicate the 
theoretical uncertainties on the default predic- 
tions. 



At NLO, the inclusive production of J/ip 
mesons with finite values of Pt has so far only 
been studied in the color-singlet channel in direct 
photon-proton [17] and proton-proton [8116] colli- 
sions and in the color-octet channel only in direct 
photon-photon collisions |14|15j . The technical 
difficulties that need to be tackled in these cal- 
culations include the treatment of UV, collinear, 
soft, and Coulomb singularities, the NRQCD op- 
erator renormalization, the mass factorization, 
and the analytic evaluation of five-point one-loop 
integrals with UV, soft, and Coulomb singular- 
ities. As for the real corrections, the phase- 
space slicing method can be employed to demar- 
cate the regions of phase space containing soft 



and collinear singularities from the hard regions, 
where the phase-space integrations can be car- 
ried out numerically. One must, of course, ver- 
ify that the combined result is, to very good ap- 
proximation, independent of the choices of the 
cut-off parameters over an extended range of val- 
ues. One usually works in dimensional regu- 
larization in connection with the MS renormal- 
ization and factorization schemes, so that the 
NLO result depends on the QCD and NRQCD 
renormalization scales /x and A, respectively, and 
on the factorization scale M connected with the 
collinear splitting of the incoming photons into 
massless qq pairs. While the y, and A depen- 
dences are formally canceled up to terms beyond 
NLO within direct photoproduction, the M de- 
pendence is only compensated by the LO cross 
section of single-resolved photoproduction. By 
the same token, the strong M dependence of the 
latter is considerably reduced. This is of crucial 
phcnomcnological importance, as may be appre- 
ciated by observing that, at LO, the overwhelm- 
ing bulk of the cross section of prompt J/ip pro- 
duction in two-photon collisions is due to single- 
resolved photoproduction. The naive K factor 
of direct photoproduction, defined as the NLO 
to LO ratio, turns out to be very substantial at 
large values of Pt because fragmentation-prone 
channels start to open up at NLO. In fact, at 
Pt ^> 2m, the Pt distribution of direct photo- 
production at NLO is rather similar to the LO 
result of single-resolved photoproduction both in 
shape and normalization. 

3. Inclusive quarkonium production in 
ultra-peripheral collisions at the LHC 

We now turn to the calculation of quarkonium 
production through photon-photon collisions in 
ultra-peripheral collisions of protons and heavy 
ions at the LHC. We focus on the production of 
quarkonia with direct photons, as only these pro- 
cesses have so far been implemented in the tree- 
level quarkonium amplitude generator MadOnia 
[27] , This Monte Carlo generator uses a sim- 
ple method to automatically evaluate arbitrary 
tree-level amplitudes involving the production or 
decay of a heavy quark pair QQ in a generic 
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2S+1 Lj' 8 ^ state, i.e. the short distance coefficients 
appearing in the NRQCD factorization formal- 
ism. The approach is based on extracting the rel- 
evant contributions from the open heavy quark- 
antiquark amplitudes through an expansion with 
respect to the quark-antiquark relative momen- 
tum and the application of suitable color and spin 
projectors. Several checks are performed to val- 
idate the calculation, including those of gauge 
invariance and numerical cancellation of contri- 
butions bound to vanish due to symmetry condi- 
tions. As a typical set for our numerical input, we 
use a J/i/j-meson mass of mj/^ = 2m c = 3 GeV, 
electromagnetic and strong coupling constants of 

a = 1/137 and a s (fj, r ) with fi r = y/ + Pj, , 
and color-singlet and color-octet operator expec- 
tation values of (O ( 3 s[ 1] )) = 1.16 GeV 3 and 

(O ( 3 s[ 8] )) = 1.06-1(T 2 GeV 3 . 

3.1. Photon flux 

A reasonable analytic approximation for the 
photon flux in ultra-peripheral collisions of pro- 
tons and/or heavy ions A and B is obtained by 
integrating the differential flux in impact param- 
eter space over radii larger than Ra + Rb, where 
the proton radius is R p = 0.6 fm and the radii 
of heavier nuclei are Ra,b ^ Ro(A,B)^ 3 with 
Ro = 1.2 fm and A, B the nucleon numbers. The 
integrated flux is then given by 

J~i/a{x) = — — 



wifo(w)ifiH (2) 



Here Z is the proton or nuclear charge, x — k/E 
is the fraction of the beam energy E of the pro- 
tons or heavy ions with mass m carried by the 
almost real photon with energy k and virtuality 



-q 2 < l/R 2 AB , if(o,i) are the th - and l st -order 



modified Bessel functions of the third kind, and 
u = xtti^Ra + Rb)- The beam energies and corre- 
sponding expected luminosities for the LHC can 
be found in Tab. 1 of [23J. Note, however, that 
for protons we do not use the spectrum above but 
rather the one in Eq. (D.7) of [28] , both for pp 
and pPb collisions. In the latter case, this gives 
us a first reasonable approximation of the photon 



flux, but the difference in the minimum value of 
the impact parameter (2i? p vs. R p + Rpb) should 
be taken into account in the future. 

3.2. Numerical results 

We present here our numerical results for 
J/V'-production in direct photon-photon colli- 
sions for three different hadron-beam combina- 
tions: proton- proton (pp) collisions at acenter-of- 
mass energy of yfs = 14 TeV, proton-lead (pPb) 
collisions at a nucleon-nucleon center-of-mass en- 
ergy of y/s = 8.8 TeV, and lead-lead (PbPb) colli- 
sions at a nucleon-nucleon center-of-mass energy 
of yfs = 5.5 TeV. The expected luminosities listed 
in Tab. 1 of [23] are 10 7 mb" 1 s" 1 , 420 mb" 1 s" 1 
and 0.42 mb" 1 s _1 . 

In order for the transverse momentum of the 
produced J/ip mesons to be balanced, it must be 
produced in association with at least one other 
particle. The direct production of a color-singlet 
with one or two gluons is forbidden by color 
conservation and charge conjugation. The same 
holds for light qq pairs, which originate from 
gluon splitting. The dominant channels arc thus 
the associated production of a color singlet with a 
photon or an additional cc pair or the associated 
production of a color octet with a gluon. Note 
that in principle these final states are experimen- 
tally distinguishable. 

The transverse-momentum spectra for J/ip 
mesons produced in pp, pPb, and PbPb collisions 
are shown in Fig. 2 (top, center, and bottom) for 
central rapidity y = 0. All cross sections have 
been multiplied by the branching ratio of 5.88 % 
for J/tp — > As one can see, the cross sec- 

tions for PbPb collisions are more than three or- 
ders of magnitude larger than the ones for pPb, 
and six orders of magnitude larger than those for 
pp. Indeed, the two-photon spectra for heavy ions 
are enhanced by factors of Z 2 for pPb and Z 4 
for PbPb up to corrections due to different min- 
imum impact parameters for an ultra-peripheral 
collisions. Note also that for instance the cross 
section for J /ipcc decreases faster as a function 
of Pt for PbPb than for pp. This can be traced 
back to the harder spectrum of 7 from p than for 
Pb. 

As we have seen above, the expected luminosi- 
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Figure 2. Transverse- momentum spectra of J/ip- 
mesons produced in association with photons 
(dotted), cc pairs (full) and gluons (dashed) in pp 
(top), pPb (center) and PbPb (bottom) collisions. 
While the first two processes produce directly a 
color-singlet, the latter produces necessarily an 
intermediate color-octet state due to color con- 
servation. 



ties for heavy-ion runs are also much smaller and 
scale, in fact, like Z 2 . Taking into account the 
run duration (10 7 sec for pp and 10 6 sec for pPb 
and PbPb), the number of events will be compa- 
rable to, if not a bit larger than in the pp runs. 
Note that ultraperipheral collisions of lead ions 
can be triggered on by requiring the detection of 
a forward neutron emitted by a lead ion, that was 
excited during the collision 23J. 

In each plot, we show three curves correspond- 
ing to the associated production of J/tp mesons 
with photons (dotted) and cc pairs (full) directly 
through an S'-wave color-singlet state and with 
gluons (dashed) through an intermediate color- 
octet S'-wave state. The associated production 
with photons and gluons receives contributions 
from exactly the same Feynman diagrams, so that 
the two curves show nearly the same shape. They 
differ only in normalization due to different cou- 
plings (electromagnetic vs. strong) and operator 
matrix elements (singlet vs. octet), the former dif- 
ference being slightly Pt dependent due to the 
running of a s that we took into account. The 
associated production of J / ip mesons with an ad- 
ditional cc pair allows, however, for different re- 
combinations of the intermediate charm quarks, 
leading to a very different and visibly much flat- 
ter Py-distribution than in the other two cases. 
In particular, this channel becomes dominant al- 
ready for P T > 2 GeV. 

4. Conclusion 

As we have seen, the mechanism of inclu- 
sive production of heavy quark-antiquark bound 
states (quarkonia) is still far from being under- 
stood. Leading-order calculations seem to indi- 
cate a large contribution from intermediate color- 
octet states in agreement with the effective field 
theory of NRQCD, but this leads to contradic- 
tions with the observed quarkonium polarization. 
Next-to-leading order corrections are known to be 
very important, but are so far only available for 
a restricted number of partonic processes. 

Proton-proton and heavy-ion collisions at the 
LHC will soon allow for studies of very high- 
energy photon collisions. We have demonstrated 
that the production of quarkonia in these colli- 
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sions will have observable rates at least in pp, 
but possibly also in pPb and PbPb scattering. 
The associated production of J ftp mesons with 
photons, gluons and cc pairs will lead to exper- 
imentally distinguishable final states, which are 
sensitive to different intermediate color states. 
Together with their discriminating transverse- 
momentum spectra, this will hopefully provide 
crucial new experimental information on this 
longstanding issue. 
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